'《Scattering of Electromagnetic Waves: Theories and Applications》第一章'

单粒子的电磁散射

本书主要研究随机分布的粒子对波的传播和散射。本章和下一章讨论单个粒子的散射。

1.基本散射参数

1.1 散射幅度与散射截面

入射波打在介电常数为的目标上,背景介电常数为。入射波沿方向,电场为:

这里,



为波数。时谐因子被忽略。

在远场区,散射场为球面波,会有衰减因子。由于麦克斯韦方程是线性的,我们可以写出:

这里垂直于散射方向,而被称为从的散射幅度。

Fig 1.1.1. 平面电磁波的散射

磁场满足:

这里为波阻抗(wave impedance)。坡印廷矢量表示单位面积的能流:

类似的,对于散射波,磁场为

坡印廷矢量:

利用(1.1.4)和(1.1.8),有

Fig 1.1.2. 微分立体角

考虑一个散射方向上的微分立体角,在球坐标系下:

在距离为处,该立体角所占的表面积为:

穿过的能量为:

将(1.1.9)代入(1.1.12):

考虑入射波的坡印廷矢量,从(1.1.6)中,我们有

(1.1.14)的量纲为面积,因此定义散射截面

对比(1.1.14)和(1.1.15),有

对(1.1.14)积分,得到

散射功率为

其中为散射截面:

散射截面和几何截面

几何截面是其投影在垂直于入射方向的平面上的区域。因此,从几何光学的视角来看,粒子“截获”的功率是几何截面和入射坡印廷矢量幅度的乘积:

我们可以比较,以及。令为物体的尺寸(物体内部两点之间的最大距离)。当远小于波长,瑞利散射理论表明:

其中表示幅度的阶数,所以

这表明粒子与波长相比很小时,粒子散射的功率远远小于

在短波的限制下,,那么

这被称为是几何光学极限。散射截面也取决于之间的差距。在弱散射体中,,我们有

这里表示正比。(1.1.24)是Born近似的结果。

吸收截面

粒子可以从入射电磁波中吸收能量,令

利用欧姆定律,吸收的功率为

这里代表粒子的内部场,,并且(1.1.26)积分是粒子体积的三维积分。吸收截面定义为

粒子总的截面为

粒子的反射率(albedo)为

因此。反射率是散射截面在总截面里面占比的一种度量。

1.2 散射幅度矩阵

我们接下来考虑极化。入射波电场垂直于传播方向。有两个垂直于的线性无关的矢量。我们称它们为

这里是正交归一矢量,遵循右手定则。

散射波类似,令形成正交系统,那么

Fig 1.1.3

散射场分量相对于线性的。该关系能够用矩阵表示:

极化描述的正交归一系统

有两个通常选择的正交归一系统:描述粒子散射。

A.基于散射平面的系统

之间的夹角为形成的平面为散射平面,令

这些时为垂直于该平面的单位矢量。并且有

根据正交性


在该1-2坐标系中,散射幅度矩阵仍然遵循:

该坐标系的优势在于对于具有对称性的粒子,散射幅度能够呈现简单的形式。缺点在于取决于散射方向。

一些有用的关系:


B.垂直和水平极化

在地球物理探测和地球遥感探测中,一般选择垂直于地球表面的垂直轴。这样我们可以得到垂直极化和水平极化,它们与构成一个正交系统。我们选择:

它垂直于,然后有

垂直极化的别名包括TM极化、平行极化和极化。水平极化的其他名称包括TE极化、垂直极化和极化。

Fig. 1.1.4 球坐标系下的入射与散射方向矢量

如果被角度确定在球坐标系里:

类似的,如果描述,那么

2.Rayleigh散射

2.1 单个小粒子的Rayleigh散射

在Rayleigh散射中,粒子尺寸远小于波长。在这种情况下,粒子内部会产生一个振荡偶极子,偶极矩为。偶极子辐射的结果就是散射场。

的远场辐射为:

令粒子内部电场由表征,其中下标表示internal。单位体积粒子的电极化矢量为:

电极化强度是指电介质在某一点单位体积内的电偶极矩。具体来说,电极化强度是一个矢量,它描述了单位体积内所有分子偶极矩的平均密度。当电介质处于外电场中时,分子内的正负电荷中心会发生相对位移,形成电偶极子,这种现象称为位移极化。电极化强度的大小与外电场强度以及介质的极化率相关,对于各向同性的线性电介质,电极化强度与外电场呈线性关系。电极化强度的散度反映了极化电荷体密度的负值,而其旋度为零,这意味着电极化强度是一个保守场。

在电磁学中,电极化强度是一个重要的物理量,它帮助我们理解电介质在外电场作用下的行为。通过电极化强度,我们可以计算出极化电荷的分布,进而分析电介质对电场的响应。电极化强度与电位移矢量密切相关,是为了解决含有电介质的静电场问题而引入的辅助场,它简化了麦克斯韦方程组的求解过程,使得在处理电介质问题时,可以忽略极化电荷的影响,直接与自由电荷关联。

电极化强度可以被定义为。而对于线性、各向同性的电介质,电极化强度满足

在Rayleigh散射中,内部场是一个常矢量。粒子的偶极矩是

其中是粒子的体积。在(1.2.1)使用(1.2.2)和(1.2.3)得到

使用(1.2.4),散射幅度矩阵可以由内部场和入射场确定。

在Rayleigh散射中,根据(1.1.26),粒子的吸收功率为

2.2 球的Rayleigh散射

考虑一个球,其半径,位于原点。粒子的内部场为:

设任意一点的总场为,该点距离球心为,高斯定理告诉我们

从而有

如果内部的球产生了电极化,那么可以分成正电球和负电球两个球,二者的球心之间有一个小距离。假定球内部有一个点,距离负电球球心距离为,距离正电球球心距离为,那么有

考虑到总电场是外电场和极化电场的叠加:

从而

考虑到线性均匀介质的关系:

因此

即有

所以对应到这里,由于粒子的介电常数为,入射电场为,所以粒子内部的总场为

内部场平行于入射场。然后利用(1.2.6)和(1.2.4),以及,有

这里

从(1.2.7)中,


我们可以用散射平面系统计算得到散射幅度矩阵。首先令,所以。那么

从(1.2.9)中,可以看到。所以

接下来,令,因此,从而。从(1.2.9)中,可以看到,。所以,从而散射幅度矩阵为

接下来,我们采用正交水平/垂直极化单位矢量。

为了得到,我们令,从而。从(1.2.9)中,我们有
因此


同理,为得到,令,有


为计算散射截面,假定入射波是水平极化的,散射功率由水平极化和垂直极化散射波在全部散射角度的总和的积分确定:

如果我们与几何截面比较的话,我们可以得到

由于,散射截面远小于几何截面。当然,该式仍然遵循(1.1.21)式,

吸收散射截面

小粒子吸收的功率为

吸收截面为

对于球的情况,将(1.2.6)和(1.2.16)代入(1.2.17)

如果我们求散射截面与吸收截面的比值,那么有

因此的相对幅度取决于

2.3 椭球的Rayleigh散射

为椭球长度的半轴,它们远小于,三个半轴形成的坐标系为。那么由入射场产生的内部电场为

其中




三个积分的和满足:

对于旋转椭球,,那么这些积分可以写成解析形式。如果为旋转轴,那么

对于长球体,即

其中

为椭球率。

对于扁球体,有

其中

对于薄圆盘,我们可以将其视为扁球体的特例,令,有



将(1.2.32)-(1.2.33)代入(1.2.20),得到

对(1.2.34)的解释是,入射电场的切向分量穿过薄圆盘,而入射电场的法向分量穿过圆盘时变为了倍。

2.4 散射并矢/散射矩阵

从(1.2.20)出发,很容易定义散射并矢:

利用(1.2.4)和(1.2.20),我们得到

所以

由于

所以




双站散射截面定义为




在后向散射时,

因此

以及

注意(1.2.43)中没有负号,而(1.2.41)和(1.2.42)中有负号。

例1

对于球体,,其中是相对介电常数。,此时

此时有

在后向散射方向:



因此



所以球体的后向散射的特征是,没有交叉极化。

如果入射波为右旋极化波,那么,此时,且。代入(1.2.47)得到,且,因此,散射波是左旋圆极化的。

例2

对于薄圆盘的情况

之后

其中

,后向散射情况。此时。那么。如果,那么。因此后向散射交叉截面的极化依赖性包含物体的形状信息。

3 散射的积分表示与Born近似

3.1 散射幅度的积分表达

本节中,我们将推导积分表达。

Fig 1.3.1 有限大小任意形状的散射体。

回想(1.2.4),单个粒子的远场散射场

考虑一个散射体,任意尺寸、任意形状、非均匀介电常数,如图1.3.1所示。令处的微分体积元,具有介电常数,此时对远场散射场的贡献

其中是源点和观测点之间的距离。如果观测点足够远,那么

我们估计(1.3.2)中的相位项:

幅度上有

对散射场进行远场近似,并对体积进行积分:

(1.3.6)中的表达式是任意大小、形状和不均匀介电常数粒子的精确远场散射振幅。然而内部场是一个未知数。为了严格计算,必须求解Maxwell方程组。例如,取决于粒子不同部分的之间的相干波作用。

除了使用内部电场,我们可以使用表面的切向电场和磁场。给定切向电场和磁场,远场辐射场为

(1.3.7)的导数可以在第2章Section2找到。

3.2 Born近似

Born近似包含用入射场近似内部场:

这适用于以下情况:

也就是说,当粒子的介电常数接近背景介质的介电常数时,粒子被称为是微弱的(tenuous)。注意Born近似并不需要假定粒子很小的假设。它假定粒子和背景介质的介电常数很小。令入射波为平面波,这样

将(1.3.8)和(1.3.10)代入到(1.3.6),得到

在(1.3.11)中,我们将积分上下限拓展至全空间,因为粒子外有,且被积函数为0。令

那么

其中

这里是粒子体积。矢量是入射波矢和散射波矢的差。

Fig. 1.3.2. 入射、散射波数矢量的关系



那么

方程(1.3.16)是三维傅里叶变换的形式。因此如果在很多的下测量,那么可以利用三维逆傅里叶变换被获取(retrieved)。该方法被用于非均匀粒子的层析成像。

然而,为了利用的三维逆傅里叶变换重构,根据(1.3.16),我们需要上的所有值。但是,由于,我们能获取的最大值。所以我们只能知道在半径为空间球形体积内的值。重建后的将会被限制在低通版本,即空间频率

如果之间的夹角,那么

最大位于,后向散射,此时;最小位于,前向散射,此时

散射场的极化为,仅仅取决于散射方向和入射极化。因此Born近似中的极化依赖性是不感兴趣的,因为它不取决于粒子的特性。

例:

对于一个半径为的球体和均匀介电常数,我们有

为了积分,不失一般性地,我们沿选择轴,此时。那么,在球坐标系下积分,有

具有处取最大值,此时。最大值为。之后随着增加,其振荡包络线逐渐减小。散射角度展宽

角度展宽意味着Born近似预测,大部分散射能量都在角度范围内。对于较大的,我们有

如果我们使用散射平面坐标系表示散射幅度矩阵,并记,然后使用(1.3.13),得到

对于球的情况,我们有

其中


由(1.3.18)给出

4 平面波,柱面波,球面波

标量波方程为

矢量电磁波方程为

本节中,我们将讨论这两个方程的解,使用的坐标系包括笛卡尔坐标系、柱坐标系、球坐标系。

4.1 笛卡尔坐标系:平面波

标量平面波解为

其中

为传播矢量。将(1.4.3)代入(1.4.1)得到

对于矢量电磁波方程,令

将(1.4.6)代入(1.4.2),得到

将上式标量点乘得到

从(1.4.7)和(1.4.8)中,我们有

从而得到(1.4.5)。

4.2 柱面波

在柱坐标系下,标量波方程为



在柱坐标系下该方程的正则解(regular solution)是

其中为整数,在(1.4.12)中,是Bessel函数,其阶数为代表正则的意思。外行柱面波为

其中是第一类Hankel函数。

函数在的行为是


且对于


函数在的情况为


对于负指标,则有


以上是对标量波方程的讨论。而对于矢量波,情况则复杂得多。矢量柱面波方程为,



当没有符号时,Bessel函数将被第一类Hankel函数取代。注意满足(1.4.2)的矢量波方程,而是无旋的,仅仅对源区场做出贡献。将与有关的项与的指数项分离便于我们后续操作。可以写为:


其中,


我们接下来试图找到平面波和柱面波之间的变换。令


我们有数学等式:

所以

利用,以及(1.4.19)

其中是水平极化矢量。因此我们用矢量柱面波表达了电磁平面波矢量:

两边取旋度,得到

然后

其中是垂直极化矢量。所以对于一个总矢量平面波,柱面波矢量为:

4.3 球面波

在球坐标系下,标量波方程为:

标量波方程实际上就是亥姆霍兹方程,它的解呈现出这样的形式(球贝塞尔方程/半整数阶贝塞尔方程):

进一步地,也可以分解为连带勒让德方程和三角函数:

其具有外行波解:

其中,。在(1.4.35)中,是第一类的球Hankel函数,且为连带Legendre函数。

这里用是因为球Hankel函数代表满足亥姆霍兹方程的球函数的分量是外行波。(是驻波,是内行波)。

我们同样可以定义正则波方程:

其中代表正则含义,为球Bessel函数。前缀“正则”用于表示代替,而时为有限值。连带Legendre函数定义为

这是罗德里格斯公式。

其对于成立,即表示(1.4.37)对的正负值均成立。正负值之间的关系满足:

前几阶连带Legendre函数如下:

对于

对于



对于





前几阶球Bessel函数为



对于球Hankel函数,我们有



由于我们仅仅使用了第一类球Hankel函数,我们将省略上标(1)。对于小段


在求复变函数的极限时,有时必须保留实部的主导项和虚部的主导项(leading term)。因此我们可以得到(1.4.44b)的小参数近似,尽管实部远小于虚部。对于的区域,

球谐函数定义为:

其具有正交关系:

完备关系为:

平面波的球波分量展开为:

其中表示描述方向的角变量。

矢量球谐函数包括三类矢量方程



注意对于从0开始,而对于从1开始。对于矢量球谐函数,其正交关系满足:

其中


三类正则矢量球面波定义为


对于三维亥姆霍兹方程,如果的解,那么任取一个无旋矢量,这三个矢量函数均满足三维亥姆霍兹方程

我们可以从平面直角坐标系下的电场去理解为什么会有三个解。取

我们这里可以看到,是沿传播方向的结果,其它的是垂直传播方向的,可以被理解为电磁波的极化。

在(1.4.55)-(1.4.57)中,


在(1.4.55)-(1.4.57)中,没有前缀的矢量球面波中的替代。我们将其记为,而对于满足方程:

满足


矢量球面波函数可以表示为矢量球谐函数的积分形式:



利用(1.4.63)-(1.4.65),可以得到

其中是与方向一致的角度变量。平面电磁波可以通过(1.4.66)和入射电场的点积,而被表示为一系列球面波。令入射波传播方向,且,那么

为了获得沿方向的传播的平面电磁波,我们令(1.4.66)中,并与极化矢量做点乘。

我们注意到:

该式除了以外都成立。因此,利用(1.4.66)

外行矢量球面波的远场解为


在Mie散射中,结果经常被表示为,它们与Legendre函数之间的关系为


时,


之后我们可以将矢量球谐函数表示为


而对于


对于,不失一般性,我们可以令,从而



5 声散射

声波是标量波的一个例子,本节将简要介绍声波传播。

运动方程为

其中是物质密度,是速度,是压强。导数是随流体运动的坐标系统的时间变化率,其通常被称为拉格朗日描述。另一方面,是一个空间中的固定点的时间导数,被称为欧拉描述。

我们有

质量守恒律有:

为了获得声方程,我们需要将(1.5.1)和(1.5.3)线性化。令



这里,我们假定为0,即流体是零阶静止的。

材料的构成关系是线性相关

常数的值取决于材料,它是声波在介质中的速度。

利用(1.5.1)式,我们得到了

将(1.5.8)式平衡(Balancing)至一阶,得到:

从(1.5.3)中,我们得到

将(1.5.10)平衡至一阶,得到:

将(1.5.7)中的本构方程代入(1.5.11),得到

从(1.5.9)和(1.5.12)中,我们可以得到的波动方程,其具有波速

两个介质之间的边界条件是:

  • 压强的连续性
  • 速度的法向分量的连续性

6 球、柱体、圆盘的散射

6.1 Mie散射

入射沿+z方向,散射体为球,介电常数为$\epsilon_p$

Mie散射是半径为,介电常数为的球体对电磁波的散射。令,且球位于原点。

由于球的对称性,很方便使用散射平面正交坐标系表示散射幅度,以得到:

我们令入射波沿方向,因此。同时,令观测方向位于平面,即,那么

入射场由(1.4.68)给出,

为了求解边值问题,我们令散射场为

并令内部场为

注意,内部场满足波数为的矢量波方程。边界条件是的连续性,其正比于。在的位置

平衡(其实就是比较两边的分量)的分量,得到


下一步,我们在处匹配

我们有,以及,将替换为即可得到内部场。因此,利用(1.6.11)和(1.6.5)-(1.6.7),有


求解(1.6.9)和(1.6.12),得到。求解(1.6.10)和(1.6.13)得到。我们发现

其中

并且

其中

我们可以充分利用朗斯基行列式(Wronskian):

其中是球Neumann函数。且


从(1.4.56)-(1.4.57)和(1.6.6)中,在远场范围内,

将(1.4.72)-(1.4.75)代入至(1.6.18),并对求和,可以得到

为了获取,我们利用(1.6.1)-(1.6.4)式,并令,从而。将这些代入(1.6.19),我们得到

因此,由于,根据(1.6.20),我们有


其中

为了获得,我们令,从而,这时(1.6.19)会变为

因此


其中

散射幅度矩阵为

为了得到散射截面,我们令

这表明。令

那么,并且



所以

其中,


然后使用(1.6.28)

所以

并且散射截面能够通过在的立体角上对总散射功率进行积分计算:

利用(1.6.23)和(1.6.27),代入(1.6.38),的正交关系满足